Magnetoelectric Coupling and Electric Control of Magnetization in 
Ferromagnet-Ferroelectric-Metal Superlattices 
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Ferromagnet-ferroelectric-metal superlattices are proposed to realize the large room-temperature 
magnetoelectric effect. Spin dependent electron screening is the fundamental mechanism at the mi- 
croscopic level. We also predict an electric control of magnetization in this structure. The naturally 
broken inversion symmetry in our tri-component structure introduces a magnetoelectric coupling 
energy of PM 2 . Such a magnetoelectric coupling effect is general in ferromagnet-ferroelectric het- 
erostructures, independent of particular chemical or physical bonding, and will play an important 
role in the field of multiferroics. 



Ferroelectricity and ferromagnetism are important in 
many technological applications and the quest for mul- 
tiferroic materials, where these two phenomena are in- 
timately coupled, is of significant technological and fun- 
damental interest [J, H, H 0, S I3| • I n general, ferroelec- 
tricity and ferromagnetism tend to be mutually exclusive 
or interact weakly with each other when they coexist in 
a single-phase material @- Increasing the spin-orbit in- 
teraction of the electrons 0] or strategically designing 
for magnetic and electric phase control 0, Q may en- 
hance the magnetoelectric (ME) coupling effect in a sin- 
gle phase multiferroic material. Practical applications of 
the ME effect, however, remain hindered by the small 
electric polarization and low Curie temperature 0, 0, 0]. 

Artificial composites of ferroic materials may enable 
the room-temperature ME effect, since both large and 
robust electric and magnetic polarizations can persist 
to room temperature. Two types of ME coupling at a 
ferromagnet(FM)-dielectric interface have been reported, 
one employing the mechanical interaction]^, [T3| or chem- 
ical bonding [Tl| and the other mediation by carriers 
(screening charges) [HI, [l3|, El, EH- The role of elec- 
trostatic screening in ferroelectric capacitors has been 
studied by macroscopic models [16J. Recently, ab-initio 
studies of nanoscale ferroelect ric( FE) capacitors [13, HI] 
and FE tunnel junctions [Til, [20, El have further con- 
firmed that electrostatic screening is the fundamental 
mechanism at the FE/normal metal (NM) interface. In 
this letter we propose a strategy of achieving robust ME 
coupling in a tri-component FM/FE/NM superlattice. 
The additional magnetization, caused by spin-dependent 
screening [lj| EH, will accumulate at each FM/FE inter- 
face. Due to the broken inversion symmetry between the 
FM/FE and NM/FE interfaces, there would be a net 
additional magnetization in each FM/FE/NM unit cell, 
unlike the symmetric structures discussed in the previous 
work. The addition of magnetization in this superlattice 
will result in a large global magnetization. 

The tri-component superlattice is illustrated in Fig.l 
(a). When the FE layer is polarized, surface charges are 
created. These bound charges are compensated by the 
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FIG. 1: (Color online). (a) Schematic illustration of 
FM/FE/NM tri-component superlattice. (b) The distribu- 
tion of charges and induced magnetization (green) calculated 
by our theoretical model. A and B are two different choices 
of the unit cell. The directions of arrows indicate the mo- 
tions of positive (red) and negative (blue) charges across the 
boundary of the unit cell A. (c) Electrostatic potential profile. 



screening charge in both FM and NM electrodes. In the 
FM metal, the screening charges are spin polarized due to 
the ferromagnetic exchange interaction. The spin depen- 
dence of screening leads to additional magnetization in 
the FM electrode as illustrated in Fig. 1(b). If the density 
of screening charges is denoted as r\ and the spin polar- 
ization of screening charges as £, we can directly express 
the induced magnetization per unit area as 



AM 



V 



Cms- 



(1) 



As this effect depends on the orientation of the electric 
polarization in FE, the magnetoelectric (ME) coupling is 
expected. 
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Before starting specific calculations, let's consider two 
simple cases. (1) In an ideal capacitor where all the sur- 
face charges reside at the metal (FM or NM)/FE inter- 
faces, the density of screening charge rj reaches its maxi- 
mum value, 77 = Po, where Pq is the spontaneous polar- 
ization of the FE. This results in a large induced magne- 
tization • C/^b)- (2) In half metals, there is only one 
type of carriers that can provide the screening. If a half 
metal is chosen to be the FM electrode, the screening 
electrons will be completely spin-polarized. In this case, 
a large induced magnetization is also expected, AM— 

Induced Magnetization from screening charges. 

For simplicity, we will first consider the case of zero bias, 
as illustrated in Fig.lc. Here, the following assumptions 
are made. (l)The difference in the work function be- 
tween FM and NM is ignored. (2)To screen the bound 
charges in FE, the charges in metal electrodes will accu- 
mulate at the FM/FE side, and there is a depletion at the 
NM/FE side. In this process, the total amount of charge 
is conserved, however, the spin density is not conserved 
because of the ferromagnetic exchange interaction in the 
FM metal. 

As shown in Fig.l, the local induced magnetization, 
defined as 6M(x) = [Sn^ (x) — 8n^(x)]fiB , is a function 
of distance from the interface x. Here, 5n cr {x) is the 
density of the induced screening charges with spin a. 
Zhang 12] considered the FM/dielectric interface within 
the linearized Thomas-Fermi model and derived two cou- 
pled equations relating the local induced magnetization 
5M(x) and screening potential Vq(x), 



by integrating SM(x) over the FM layer, 



SM(x) = -M^-eV (x) 



d 2 Va(x) 
dx 2 



-V Q (x) 



(2) 



The screening length in the FM electrode is defined 



-1/2 



as X FM = j , where N = 

N^+N l is the total density of states, M = (iV T - N l )u B 
can be thought of as the spontaneous magnetization, Eq 
is the vacuum dielectric constant and J is the strength 
of the ferromagnetic exchange coupling in the FM layer. 

We solve the above equations for our unit cell and ob- 
tain 



V (x) 



r)*NM p-{x-t FE )/\ffu 



-t FM < x < 0; 

tpE < X < tpE + tpjM- 

(3) 

where 77 is the density of screening charges, Xfm(nm) 
is the screening length of FM(NM) electrode, and tpM, 
tpE and t n m are the thickness of FM, FE and NM layer, 
respectively. From above equations, we see that the local 
induced magnetization SM(x) decays exponentially away 
from the FM/FE interface. This distribution is identical 
to that of screening charges, because in our model the 
effective interaction J in FM is assumed to be a constant. 
The total induced magnetization AM can be calculated 



AM = / 5M(x) 

'FM layer 



N + JN* - J(M /ti B y 



(4) 

The effective spin polarization of screening electrons can 
then be written as 



c = - 



M 



N + JN$-J{M /u B y (5) 

We have considered the induced magnetization in 
FM/FE/NM tri-component superlattice with several FM 
electrodes, i.e., Fe, Co, Ni and CrC>2. Detailed parame- 
ters and calculated values of AM are listed in Table I. 
The magnitude of AM is found to depend strongly on 
the choice of the FM and FE. Among the normal FM 
metals (Ni, Co and Fe), the largest AM is observed in 
Ni for its smallest J and highest spontaneous spin polar- 
ization Mo/ (ibNq. On the other hand, we also predict a 
large AM for the 100% spontaneous spin polarization in 
half- metallic Cr02. 

To confirm the validity of our model, we perform first- 
principles calculations of the Fe/FE/Pt superlattice [23| . 
The calculations are within the local-density approxi- 
mation to density-functional theory and are carried out 
with VASP We choose BaTi0 3 (BTO) and PbTi0 3 
(PTO) for the FE layer. Starting from the ferroelectric 
P4mm phase of BTO and PTO with polarization point- 
ing along the superlattice stacking direction, we perform 
structural optimization of the multilayer structures by 
minimizing their total energies. The in-plane lattice con- 
stants are fixed to those of the tetragonal phase of bulk 
FEs. Fig. 2 shows the calculated induced magnetic mo- 
ment relative to that of bulk Fe near the Fe/BTO in- 
terface when the polarization in BaTiOa points toward 
the Fe/BTO interface. It is evident that the induced 
moments decay exponentially as the distance from the 
interface increases. This result is in line with our model 
for the magnetization accumulation in the FM at the 
FM/FE interface. A numerical fitting of the exponen- 
tial function yields a screening length of ~ 0.7 A for 
the Fe/BaTiOa/Pt structure. This value is comparable 
to the screening length parameters calculated using the 
theoretical model as shown in Table I. 

Using our theoretical model we also calculate the ME 
coupling coefficient r, which is defined as the ratio of the 
magnetization change 2(1qAM/ Xfm to the coercive field 
Vc/(l*tpE), where Vc is the coercive voltage and I is the 
number of the unit cell. The r values listed in Table I 
reach as large as 0.015 G cm/V in Ni. For comparison, r 
of about 0.01 G cm/V arising from the chemical bonding 
between Fe and Ti atoms is predicted for the Fe/BaTi03 
bilayer [ll[. The ME coefficient measured in epitaxial 
BiFe03-CoFe204 columnar nanostructures [10( is also of 
0.01 G cm/V. We should point out that in our calcula- 
tion the coercive field is assumed to be 200 kV/cm, if we 
choose the coercive field of 10 kV/cm same as Ref. [TJ, r 
will be 20 times larger than those listed in Table I. There- 
fore, the ME effect arising from spin-dependent electron 
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TABLE I: Parameters extracted from band structures of Ni, Co, Fe from Ref . [l2l] . Cr02 from Ref.[%|. Here, AM is the value 
at V a ~Vc, where V a is the applied bias and Vc is the coercive bias. 



FM 


J (eV-nm 3 ) 


No (eV^-nm^ 3 ) 


Mq/ubNo 


Xfm (A) 


AM (/j s -nnr 2 ) 


r (G cm/V) 


Ni 


0.65 


1.74 


-79.3% 


0.9 


-0.280 


0.015 


Co 


1.25 


0.89 


-58.4% 


1.5 


-0.126 


0.004 


Fe 


2.40 


1.11 


56.8% 


1.3 


0.078 


0.003 


CrQ 2 


1.8 


0.69 


100% 


1.7 
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FIG. 2: (Color online). Layer-resolved induced magnetic mo- 
ment of Fe near the interface between Fe and BaTiOs in the 
Fe/BaTiOa/Pt superlattice. Solid line is the fitted exponen- 
tial function for the induced moment as a function of the 
distance from the interface. 



screening in FM/FE/NM tri-component superlattice can 
be much larger than in other composite multiferroics. 

What is the source of this large ME effect? In fact, the 
magnetoelectric effect discussed in this letter is not the 
usual bulk magneto-electric coupling at all. Spontaneous 
electric polarization in FE results in the induced surface 
charge. In turn, this produces the screening charges of 
density 77. These screening charges are polarized with the 
polarization (. Therefore, it is the amount of screening 
charges and polarization that determine the magnitude 
of the ME effect (AM and r). If we expand the in- 
duced magnetization AM in Eq.(4) as power series in 
order parameters Pq and Mq (spontaneous polarization 
and magnetization), we obtain 

AM oc -PcMo + high-order terms in Po and Mo (6) 

The higher order terms in Eq.(3) vanish exactly in the 
following limiting case: the screening length Xfm ~ y 
and spin polarization Q — > ±100%. In general, the lead- 
ing term in Eq.(6) is linear, which is consistent with the 
computational result of Ref. [l3[ . It is also clear that this 
effect depends on the magnetization of the ferromagnetic 
metal. 

First-principle calculations confirm the central conclu- 
sion that the ME coupling in the tri-component system 
is linear in polarization of FE. We compare a superlattice 



with BTO and that with PTO. The induced magnetiza- 
tion difference is 3.6 times larger in the the case of PTO. 
This ratio is almost exactly that of the bulk spontaneous 
polarization of BTO and PTO. 

Electric control of magnetization. So far, we have 
discussed the magnetoelectric coupling effect in the case 
of no external bias. A natural question is what happens 
to AM when external bias V a is applied. In this case, 
the electric polarization P will have two parts: the spon- 
taneous polarization Po and induced polarization. The 
equation determining P is obtained by minimizing the 
Free energy. From the continuity of the normal compo- 
nent of the electric displacement, we find equation relat- 



ing 77 and P: rj = ( 



P-u 



)/0 



f+Aj> 



EFE t n v £0 SB 

Efe is the dielectric constant of the FE layer. 



, Here, 
These 



two equations need to be solved self-consistently. The 
value of rj at a given bias can then be calculated and the 
induced magnetization AM is given by Eq.(4). 

The free energy density F includes contributions from 
the FE layer, FM layer and FM/FE interface and takes 
the form 



F = 



t FE ■ F{P) + t FM ■ F(M) + Fr(P, M) 

tFE + tFM + tjfM 



(7) 



M is the magnetization of the bulk ferromagnet, and 
here M=M because of zero external magnetic field. The 
interface energy Fj (P, M) is the sum of the electrostatic 
energy and magnetic exchange energy of the screening 
charges 



Fi(P,M) 



(Xfm + Ajvm) 2 

T) 

2e ' 



J 
2/4 



(M + AM) ■ AM. 

(8) 

For FE, the free energy density F(P) can be expressed 
as F(P) = F P + a P P 2 + (3 P P 4 + E B dP, where F P is 
the free energy density in the unpolarized state, ap and 
[3p are the usual Landau parameters of bulk ferroelectric. 
Eb is the depolarization field in the FE film. Similarly, 
F(M) can be expanded as a series in order parameter 
M, i.e., F(M) = F M + HmM 2 + v M M 4 , where F M is the 
free energy density of bulk ferromagnet, hm and vm are 
the Landau parameters of bulk ferromagnet. The calcu- 
lated induced magnetization as function of applied bias 
is shown in Fig. 3. Clearly, the electrically controllable 
magnetization reversal is realized. 

To discuss the macroscopic properties of the electric 
control of magnetization, we analyze the magnetoelectric 
coupling energy in our tri-component superlattice. For 
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FIG. 3: (Color online) AM versus V a /l for different ferro- 
magnetic metal electrodes. V a is the applied bias and / is the 
number of the unit cell. Here, the thickness of FE layer is 3 
nm. However, a thicker FE layer can be used to avoid the 
possible electron tunneling effect. 

the macroscopic average polarization to be represented 
by the electric polarization obtained for a unit cell, this 
cell needs to be chosen with special care [25j. Therefore, 
in the following calculation of total free energy, unit cell 
B in Fig. 1(b) is chosen, and 

p _ P ■ tf\E + VjtFM + tNM) /g\ 
t FE + t FM + t NM 

The macroscopic average magnetization M is 
M — ^™ + AAf 

tFE + tpM + tNM 

Considering the lowest order term of the magnetoelec- 



tric coupling, P and M can be expanded as P — c p P + 
c' p PM 2 ;M = c m M + c' m PM. Therefore, the total free 
energy (Eq.(7)) can be expressed as the power series of 
P and M, F(P, M) = F + aP 2 + f3P 4 + t iM 2 + vM A + 
\PM 2 + ■ ■ ■ . We would like to point out that biquadratic 
ME coupling P 2 M 2 is easily achievable, but is usually 
weak and is not electrically controllable. However, be- 
cause of the naturally broken inversion symmetry, the 
large ME coupling PM 2 is possible in our tri-component 
structure. 

The ME coupling in FM/FE/NM superlattices may be 
observed experimentally and may have practical applica- 
tions. Though the net additional magnetization of each 
FM/FE/NM unit cell is small, stacking several of them 
in a superlattice will result in a large overall magneti- 
zation. From Eq.(10), with a thinner metallic electrode, 
the ME effect will be larger, as long as the thickness of 
metallic electrodes are larger than the screening length, 
which is easy to achieve. 

To summarize, expanding upon the theory of spin- 
dependent screening we develop a theory of addi- 
tional magnetization in tri-component superlattice. We 
show that the additional magnetization can be electri- 
cally controlled, and is linear in FE polarization. The 
latter can be switched if the coercive voltage of the ferro- 
electric is reached. We demonstrate that an asymmetric 
FM/FE/NM structure has practical advantages over pre- 
viously discussed symmetric structure. 
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